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ABSTRACT 

We have carried out a sensitive high-resolution imaging survey of stars in the young (6-8 Myr), nearby (97 pc) 
compact cluster around 77 Chamaeleontis to search for stellar and sub-stellar companions. Given its youth and 
proximity, any sub-stellar companions are expected to be luminous, especially in the near infrared, and thus 
easier to detect next to their parent stars. Here, we present VLT/NACO adaptive optics imaging with compan- 
ion detection limits for 17 r^Cha cluster members, and follow-up VLT/ISAAC near-infrared spectroscopy for 
companion candidates. The widest binary detected is ^0"2, corresponding to the projected separation 20 AU, 
despite our survey being sensitive down to sub-stellar companions outside 0."3, and planetary mass objects 
outside 0."5. This implies that the stellar companion probability outside 0"3 and the brown dwarf compan- 
ion probability outside 0."5 are less than 0.16 with 95 % confidence. We compare the wide binary frequency 
of ?7Cha to that of the similarly aged TW Hydrae association, and estimate the statistical likelihood that the 
wide binary probability is equal in both groups to be < 2 x 10~ 4 . Even though the ryCha cluster is relatively 
dense, stellar encounters in its present configuration cannot account for the relative deficit of wide binaries. We 
thus conclude that the difference in wide binary probability in these two groups provides strong evidence for 
multiplicity properties being dependent on environment. In two appendices we derive the projected separation 
probability distribution for binaries, used to constrain physical separations from observed projected separations, 
and summarize statistical tools useful for multiplicity studies. 

Subject headings: binaries: general — stars: low-mass, brown dwarfs — stars: pre-main sequence — stars: 
planetary systems 



1. INTRODUCTION 

Mos t solar-type sta rs reside in binaries 
( Duau eimov & MavorJ Q291), yet their formation is not 
well understood. It is generally accepted that most stars 
form in groups rather than in isolation ( Adams & Mvers 
120011) . Intriguingly, the multiplicity of young stars in many 
regions seems to be systematically high er than that of 
their main-sequence counterparts (e.g. Mathieu et al. 2000, 
and references therein). The reason for this multiplicity 
overabundance is not entirely understood, although two 
principal scenarios have been proposed. One is that the 
multiplicity fraction is sensitive to the initial conditions of 
the star formation process. That would imply that different 
regions show a variety of multiplicity properties, and in- 
deed, there seems to be some evidence for the multiplicity 
fraction being anticorrelated with the stellar density of the 
region ( Patience & Duchene 2001). Another possibility is 
that young stars start out with a high fraction of multiples 
that subsequen tly are disrupted due to dynamical evolution 
(e.g. iReipurthl [2000). Star-forming regions are generally 
too dispersed for binaries to be disrupted by interactions 
between members. The alternative is that many stars form not 
just in binaries, but in unstable higher-order multiples that 
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get disrupted with time , as in deed suggested by numerical 
simulations ( Ba te et all 120021 iDelg ado-Donate et al. 2004). 
For a recent review on the current status of the field, see 
iDuchene etail J2006I) . 

The nearby (97 pc), young C6-8Mvr. Ijilinski et alJ^005t 
and references therein) 77 Chamaeleontis cluster was found by 
the X-rays emitted from its members, as revealed by ROSAT, 
together with th eir common sp ace motion, as revealed by Hip- 
parcos (Mamai ek et alJll999l) . The cluster is mainly popu- 
lated by late-type stars ( K3-M6, see Ta ble [D, and shows no 
evidence for extinction (Mamaiek et al. 2000). These prop- 
erties all contribute to make the cluster an excellent labo- 
ratory for investigating brown d warf (BP) and p l anet for- 
matio n, and their evolution (e.g. ISong et al.ll2004t IL uhman 
[2004 ]lHaisch et alJ2005ULvo et alJ2006T1jav awardha na et alJ 
2006). In particular, using adaptive optics (AO) systems on 
large telescopes, it is possible to reach high contrast-ratio sen- 
sitivities close in, to detect any wide (>0"5) companion down 
to planetary masses. 

In a previous high-resolution survey using speckle 
interfero metry and AO on 2.2-3.6 m telescopes, 
iKohler & Petr-Gotzensl J2002I) searched 13 members for 
companions (r^Cha 1-12 and 15, see Table They found 
two resolved binaries, ^Cha 1 and r^Cha 9, and one sus- 
pected unresolved binary, r^Cha 12. ??Cha 12 has also been 
suspected t o be binary due to its elevation in color -magnitude 
diagrams (Laws on et alJ200 It ILuhman & Steeghs 2004). 
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FIG. 1. — The close 0"2 binary »7Cha 1 observed with a semi-t ransp arent 
coronographic mask, and a likely background object 9" away (§ 14.11 . The 
orientation and scale are shown by the two 1" perpendicular axes.. 




FIG. 2. — Close-ups of the tight binaries r\ Cha 9 (left) and r\ Cha 12 (right). 
The box sizes are 1" on a side. (yCha 9 is a well resolved 0"2 binary, 
while 77 Cha 12 is unresolved. Fro m th e asymmetric profile, we infer that 
the rjCha 12 separation is ~0"04 (S I3.1I . North is up, and east is to the left. 

In the present paper we report on a deep AO search for faint 
companions to 17 confirmed members of the 77 Cha cluster. 
§|2]details the observations. §[3]describes analysis and results, 
including how astrometry and photometry were performed, 
how the Strehl ratio was measured and the contrast sensitiv- 
ity estimated. In §|4] we discuss the companion candidates 
and constrain the orbit of the previously observed 0."2 binary 
77 Cha 1 (Fig. [l). We then go on to determine limits on the 
companion probability, use this to discuss the apparent deficit 
of binaries in the 77 Cha cluster, and estimate the likelihood of 
binaries being disrupted by stellar encounters. §[5]contains an 
enumerated list of conclusions. 

In the Appendices, we first derive the (time-averaged) prob- 
ability distribution for the projected separation of a binary, 
given its semi-major axis. We then continue with some useful 
properties of multiplicity statistics, in particular how to es- 
timate binomial confidence intervals, and how to test if two 
outcomes are derived from the same binomial distribution. 

2. OBSERVATIONS 

We have observed 17 confirmed members of the 77 Cha 
cluster with high-resolution AO imaging and follow-up spec- 
troscopy of companion candidates. Table summarizes the 
observed target parameters and assigns a running target iden- 
tification number used in this paper. Imaging data were ob- 
tained in service mode with the AO system NACO on the 
8.2 m VLT Yepun at ESO (Cerro Paranal, Chile), during two 
different semesters (2002/2003 and 2004). Spectroscopic 
follow-up observations were obtained in service mode during 
2004 using NACO and ISAAC (on VLT Antu). 

For imaging with NACO, we used the high-resolution lens, 
giving a pixel scale of 13.26 mas pix" 1 and the field of view 
1 3 ."6 x 1 3 "6. The full-width half maximums (FWHMs) of the 
diffraction limited point-spread functions (PSFs) of the setup 
in our used photometric bands J (1.27 /im), H (1.66 /im), 
and K s (2.18/im) are 32 mas, 42 mas, and 56 mas, respec- 
tively. In the first semester, we obtained an extra exposure 
with a 5 mag neutral density filter for the stars we expected 



would saturate the array during the minimum exposure time 
of 0.3454 s. In the second semester, we instead made use of 
the new semi-transparent coronographic mask for the bright- 
est targets, that reduces the light within 0"35 radius. We mea- 
sured the contrast between the inside and outside of the mask 
to be 6.0 ±0.1 mag in H and 6.3 ±0.1 mag in K s , by observ- 
ing a binary with and without the coronographic plate. This is 
consistent with results found by the NACO instrument team 
(N. Ageorges, private communication). We used either the vi- 
sual wave front sensor (WFS) mode VIS, or the infrared WFS 
modes N20C80 or N90C10, that direct 20% or 90% of the IR 
light to the WFS, and the rest to the science camera CON- 
ICA. The zenith seeing in the V-band was better than 0"6 in 
general, though the targets, due to their low declination of - 
78°, were observed at the relatively high air mass of 1.5-2.0. 
Three examples of obtained images are show in Figs.^&El 

For spectroscopy with NACO, we used Grism 3 with the 
S27 lens, the N90C10 WFS, and a 0"086x40" slit that 
produces an //-band spectrum from 1.44/im to 1.72 /jm at 
spectral resolution R ~ 1500 and pixel scale 27 mas pix -1 x 
3.4 A pix" 1 . 

For spectroscopy with ISAAC, we used the SWS1-LR 
mode with the SH filter, and a l"x 120" slit that produces an 
//-band spectrum from 1.48 /im to 1.80 71m at spectral resolu- 
tion R ~ 500 and pixel scale 0."148 mas pix" 1 x 4.1 A pix" 1 . 

The basic reduction was done in a standard way, making use 
of the reduction pipeline in the case of ISAAC. The sky was 
estimated from the jittered observations and subtracted from 
all frames, which were subsequently corrected for cosmic rays 
and flat fielded. Since the coronographic mask is in a fixed po- 
sition on the array, half of the integration time was spent with 
the source chopped out of the field. For the spectroscopic ob- 
servations, the slit was put over both primary and companion 
candidate, and then jittered along the slit. To decompose the 
spectra of two stars on the slit, we extracted the spectra by 
fitting two-component Moffat functions (Moffat 1969) to the 
spatial profiles. The wavelength calibration made use of Ar 
lamp spectra obtained during daytime. To correct for atmo- 
spheric absorption lines, an early type (B2-B5) telluric stan- 
dard star was observed each night, at a similar airmass (within 
0.2) to the target. The science spectra were then divided by 
the telluric standard spectra, multiplied by standard star mod- 
els, and normalized. Unfortunately, the telluric line correction 
proved to not be entirely reliable, most probably due to the 
high airmass (~2) at which the observations were made. 

The use of AO and a narrow slit makes the calibration of 
NACO spectra difficult. For example, the PSF, and therefore 
slit loss, varies with wavelength. In the NACO //-band spectra 
of the resolved 77 Cha 9 binary we noticed the spectral shape 
to be somewhat steeper than the ISAAC //-band spectrum we 
have for the unresolved binary. Since the observations were 
made at high airmass, we suspect this difference might be due 
to additional wavelength-dependent slit losses, caused by at- 
mospheric differential refraction. To test this hypothesis we 
computed the expected atmospheric dispe rsion using the re- 
fraction index from Peck & Reederj ill 9721). and t he standard 
dispersion equation (e.g., equation 3 in Roe 2002), giving the 
expected dispersion of 50 mas between 1.5 /im and 1.7 /zm, 
which is a fair fraction of the 86 mas slitwidth. Projected on 
the slit orientation, the computed offset increases from 25 mas 
to 34 mas during the 10 exposures of rj Cha 9, while at the 
same time the observed spectral slope gets steeper by ~15%, 
indicating that differential refraction may indeed be signifi- 
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TABLE 1 

Targets observed in the rj Chamaeleontis cluster 



»7Cha a 


Name 


«J2000.0 


<5j2000.0 


J 


H 




SpT 


1 


EG Cha 




-78:56:45.5 


o ice I nmn 

S.155 ±o.oiy 


/.4yo ± o.04y 


/.33o ± 0.021 


K4 


z 


?7 Cha 


no ■ A 1 ■ 1 n a q 
0o.41.iy.4o 


1Q ■ C 1 ■ A Q 1 

- /o.J / .4o. 1 


5. oos ± o.oiy 


5. /zl ± 0.040 


5. /lo ± O.Olo 


B8 


3 


EH Cha 


no . a 1 .m no 

08:41:3 /.03 


- /9:03:3U.4 


1 A t /in I n (Til t 

10. 349 ± 0.023 


n c A i I n 

9.647 ± 0.022 


C\ A 1 c I nmn 

9.415 ± 0.019 


M3 


4 


EI Cha 


08:42:23.72 


-79:04:03.0 


n c o c 1 r\ a 

9.535 ±0.024 


o Tin i n n/"" i 

8.779 ± 0.061 


8.615 ±0.019 


K7 


5 


EK Cha 


08:42:2 7.11 


TO . CI. ^ -7 CI 

-/8:5/:4/.9 


10.777 ± 0.023 


10.099 ± 0.021 


9.855 ± 0.021 


M5 


6 


EL Cha 


08:42:38.80 


-78:54:42.8 


10.232 ±0.027 


9.584 ±0.023 


9.290 ±0.021 


M2 


7 


EM Cha 


08:43:07.24 


-79:04:52.5 


8.420 ± 0.024 


7.758 ±0.034 


7.635 ±0.033 


K3 


8 


RSCha 


08:43:12.23 


-79:04:12.3 


5.994 ±0.030 


5.877 ±0.038 


5.852 ±0.034 


A7 


9 


EN Cha 


08:44:16.38 


-78:59:08.1 


10.260 ± 0.026 


9.668 ±0.026 


9.335 ±0.024 


M4 


to 


EOCha 


08:44:31.88 


-78:46:31.2 


9.653 ±0.023 


8.919 ±0.063 


8.732 ±0.021 


K7 


11 


EPCha 


08:47:01.66 


-78:59:34.5 


8.729 ± 0.020 


8.025 ± 0.055 


7.655 ±0.038 


K4 


12 


EQCha 


08:47:56.77 


-78:54:53.2 


9.323 ± 0.024 


8.683 ±0.082 


8.410 ±0.031 


Ml 


13 


HD 75505 


08:41:44.72 


-79:02:53.1 


7.059 ± 0.026 


6.987 ±0.036 


6.928 ±0.023 


A1V 


15 


ECHAJ0843.3-7905 


08:43:18.58 


-79:05:18.2 


10.505 ±0.026 


9.834 ±0.021 


9.431 ±0.023 


M4 


16 


ECHAJ0844.2-7833 


08:44:09.15 


-78:33:45.7 


12.505 ±0.024 


11.976 ±0.022 


11.618 ±0.024 


M5 


17 


ECHAJ0838.9-7916 


08:38:51.50 


-79:16:13.7 


11.275 ±0.023 


10.721 ±0.022 


10.428 ±0.023 


M5 


18 


ECHAJ0836.2-7908 


08:36:10.73 


-79:08:18.4 


11.849 ±0.024 


11.277 ±0.026 


10.945 ±0.021 


M5.5 



NOTE. — Coordinates and IR magnitudes are from the 2MASS All Sky Data Release. The spectral types are from Mamaiek et al. ' 1999) for 
l-12. IHouk & Cowle\lH97l for 13. lLawson etaT]|200l for 15, and lSong et all 12001 for 16-18. 

\ 120041 . Their number 



a Numbers 1-12 coincide with the RECX numbers introduced by Mamaiek et al. 1999 ), and 1-18 
14 ("USNO Anon 1") was not observed. 



cant. The telluric standard was observed with the slit at the 
parallactic angle, and was thus not affected by differential re- 
fraction. These wavelength-dependent slit losses are of no 
consequence for the present observations, since we are inter- 
ested in the spectral difference between the two components, 
but they might be of importance for future observations using 
VLT/NACO at high airmass. 

Table|2]presents the observing log, with epochs, instrumen- 
tation, exposure time, and total integration time on source for 
each target. 

3. ANALYSIS AND RESULTS 
3.1. Astrometry and photometry of sources in the field 

For observations where multiple sources were well sepa- 
rated in the field, we measured the pixel positions on the ar- 
ray by making use of the iraf routine imexamine. To de- 
termine the relative precision of the astrometry we measured 
the separation between the tight, but well resolved, binary 
77 Cha 1 for 10 different consecutive frames. The standard 
deviation was found to be about 0. 12 pixels, or er sra =1.6 mas. 
In addition, there is a systematic image scale uncertainty of 
Csci ~ 0.01 maspix" 1 , so the estimated total separation error 
er sep depends on the separation s as <7 sep (.s) = [(<T S cis) 2 +cr^ td ] 1/ ' 2 . 
For some of the fainter sources, cr st d would be larger by a fac- 
tor of several, due to poor centering. 

The relative position angle error <tpa was computed from 
Cstd by scaling the error with separation, i.e. <Jpa(«) = 
arctan(cT st dS _1 ) ~ o- stc \s~ l rad. The position angle error is gen- 
erally dominated by the systematic error of ~ 1°, due to the 
uncertain orientation of the array. The astrometry is summa- 
rized in Table|3]&El 

Photometric measurements of AO data are complicated by 
the spatially varying PSF To estimate the level of anisopla- 
natism, we measured the PSF for the targets with multiple 
sources in the field. The FWHM of the PSF core (and thus 
the Strehl ratio) was found to be strongly dependent on the 
distance to the primary WFS. Out to a few arcseconds, the 
variation is less than 10%, but at larger separations (>7") the 
Strehl ratio can decrease by as much as a factor of 10. Part of 
the reason for this strong anisoplanatism may be that the stars 
were all observed at relatively high airmass (77 Cha never rises 



below an airmass of 1.5 from Cerro Paranal). 

Because of the dependence of the PSF on separation, 
sources separated >4" from the primary have poorly con- 
strained photometry. For those well-separated sources, we 
used a large aperture of 100 pixel radius for both the pri- 
mary and secondary, if it was bright enough; otherwise we 
fit a Gaussian to the PSF core and computed the integrated 
flux under the Gaussian. Since the Strehl ratio is so low for 
the sources at large separation from the WFS star, their PSFs 
are reasonably well approximated by Gaussians. By using the 
Gauss-fitting procedure also on the bright wide companions, 
we found the two methods to be consistent within 0.3 mag, 
which is the quoted error. As a second consistency check we 
used the relatively bright source 7" from 77 Cha 4, catalogued 
as H = 13.81 ±0.32 by 2MASS, in agreement with our esti- 
mated H= 13.9 ±0.3. 

For the close-in companions, where the PSF is roughly con- 
stant, we used aperture photometry with a radius of 4 pixels, 
and subtracted the background halo from the primary by mea- 
suring its brightness at the same separation but opposite posi- 
tion angle. The photometry is much more robust in this case. 
The error of the measured flux ratios is estimated by finding 
the consistency of multiple (consecutive) observations of the 
same target. 

As the binary 77 Cha 12 is unresolved (Fig. it requires 
special attention. We assume that the PSF core of the obser- 
vation is circularly symmetric, as is the case for the other ob- 
served central PSFs, and that the elongation to the north east 
is the result of an equal mass binary. The position angle is 
then estimated from the position angle of an elliptical Gaus- 
sian fit, and the separation from the FWHM of the fit along 
and across the major axis, by computing their difference. 

3.2. Strehl ratios 

The Strehl ratio S is defined as the ratio between the ob- 
served peak intensity P b s of the PSF, and the theoretical 
peak intensity Pjdeai of an idealized telescope (no distortion 
of wavefront, no obstructions) of the same aperture observing 
the same star, S = fobs/fldeai- Instead of computing Pjdeai for 
each observation, we estimated the ratio Rid ed \ = AdeaiMideab 
where A;d ea i is the flux (integrated intensity), that only de- 
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TABLE 2 
Observation log 









DIT 


— 

Total 


rj Cha d 


UT Date 


Instrument 


(s) 


(min) 


I 


9 004 -04 -fn 


VT<\ N rnr 

V LO 1 1 \j\JL, 




19 


2 


9004-04-04 


VIS H rnr 

V L\J 1 1 . 


3 


15 




9004-04-0? 


NQOriO 7 


20 


13 


a 

j 


900? 10 90 


VIS // 






3 


9009-10-91 


VIS // 


0.6 


14 


3 


2004-04-02 


N90C 1 K 


20 


27 


4 B 


2004-04-22 


ISAAC H sner 


200 


40 


4 


9009 11 17 

ZUUZ- 11-1/ 


VTS N wn 
v iij n in u 


i .j 


7 


4 


9009 11 17 

ZUUZ 11-1/ 


VIS // 


34 


1 3 


g 


900^-01-17 


N90CR0 hi 


5 


20 




2003-01 -21 

Z_V 'V./- ' V 7 1 1 


N20C80 H 


0.5 


14 


5 


900^-09-17 


N90CR0 N 


5 


6 




900^-09- 1Q 


N90PS0 hi 


5 


1 3 




9003-01-1 S 

V 1 1 o 


v u 1 1 x > Ly 


5 


1 .6 


5 


2001-01 -IS 


VIS // 


0.35 


21 


7 


9004-04-04 


V Ihj 11 LiUli 


20 


20 


^ 


9004 04 04 


V 1l3 /7 L U1. 




14 


9 


9004-04-09 


NQ0C10 i 


X 


13 


9 


9003-01-91 


v u 1 1 iysLs 




2 


9 


9003-01 -91 

Z.UUJ U 1 i. 1 


VIS H 


35 


1 3 


9 


9004-04-09 


NQ0C10 K 


16 


27 


9 Aab 


9004 04 03 


NTQ0C1 f/ snpr 

i>i7ULiu /7 spec 


1 20 


40 


9B 


9004-04-99 




200 


47 


10 


9003-01 -99 


V IvJ / i I'll-/ 




5 


10 


9003-01 -99 


VIS H 


35 


1 3 


10B 


9004-04-94 


TSA AC 1 /7 crwr 


200 


60 


1 1 


2001-01 -22 


VIS H ND 


1 .2 


0.4 




9003-01 -99 


VIS // 


35 


1 3 


12 


2003-01-22 


VIS // ND 


1.5 


0.5 


12 


2003-01-22 


VIS// 


0.35 


13 


13 


2004-04-03 


VIS H cor. 


10 


12 


15 


2004-04-02 


N90C10J 


16 


27 


15 


2003-01-21 


VIS H ND 


6 


4 


15 


2003-01-21 


VIS// 


0.35 


11 


15 


2003-01-22 


VIS// 


0.35 


11 


15 


2004-04-02 


N90C10 K s 


16 


26 


15C 


2004-04-22 


ISAAC // spec. 


200 


40 


16 


2004-04-02 


N90C10 H 


80 


36 


17 


2004-04-03 


N90C10 H 


10 


33 


18 


2004-04-02 


N90C10 H 


30 


33 



a Letters refer to the companion candidate observed in addition to 
the primary. 

b For NACO observations, VIS, N20C80 and N90C10 correspond 
to the used WFS, J, H, and K s are the used filters, ND means neutral 
density filter, cor. coronographic observations, and spec, spectroscopy. 
All ISAAC observation were made in the same mode (see ^2l . 



TABLE 3 
Confirmed companions 



ryCha 


Sep (") 


PA a (°) 


Band 


FR b 


1 


0.192 ±0.002 


354.0 ±0.3 


H 


1.15 ±0.02 


9 


0.204 ± 0.003 


195.7 ±0.5 


.1 


1.03 ±0.04 


9 


0.207 ± 0.003 


195.8 ±0.5 


H 


0.96 ± 0.04 


9 


0.209 ± 0.003 


196.4 ±0.5 


K- 


1.05 ±0.04 


12 c 


0.040 ±0.010 


28±4 


H 


1.0± 0.1 



a The position angle is measured from north to east. There is an additional 
systematic error of ~ 1 ° . 
b The flux ratio A/B. 
c 77 Cha 12 is unresolved, see s l3.ll 



pends on telescope aperture and wavelength. For an 8.2 m 
telescope we computed Rj = 772 arcsec" 2 , Rh = 452 arcsec" 2 , 
and Rk s = 263 arcsec" 2 , for J, H, and K s , respectively. We 
then estimated similar ratios R \, s = P b s /A ^ s for the observa- 
tions, and determined the Strehl ratio as S = /? bs/^ideai- The 
peak intensity of the observed profile was found by fitting a 
circular Gaussian function to the PSF core, and the integrated 



TABLE 4 
Likely chance alignment objects 



t) Cha 


O /// \ 

Sep( ) 


PA ( ) 


Band 


Mag 


l c 


8.61O±0.O18 


13.8 ±0.2 


H 


11.5 ± 1.0 


2 


7.997 ± 0.020 


263.6 ±0.2 


H 


19.8 ±0.5 


3 


2.084 ± 0.005 


106.5 ±0.3 


.1 


17.0 ±0.5 


3 


2.034 ± 0.005 


106.6 ±0.3 


H 


16.3 ±0.2 


3 


2.082 ± 0.005 


106.1 ±0.3 


Ks 


15.9 ±0.2 


¥ 


7.352 ±0.018 


273.3 ±0.2 


H 


13.9 ±0.3 


7 


5.611 ±0.016 


344.2 ±0.2 


H 


15.9 ± 1.0 


8 


9.443 ±0.023 


30.8 ±0.2 


H 


17.4± 1.0 


9 e 


3.582 ±0.010 


141.7 ±0.3 


.1 


15.6 ±0.2 


9 e 


3.526 ±0.010 


142.2 ±0.3 


H 


15.3 ±0.1 


9 e 


3.568 ±0.010 


141.8 ±0.3 


K s 


15.1 ±0.1 


10 


9.895 ± 0.026 


62.5 ±0.2 


H 


16.6 ±0.5 


I5 f 


2.707 ± 0.040 


73.4 ±0.8 


.1 


18.0±0.5 


15 


2.726 ±0.008 


72.1 ±0.3 


H 


17.5 ±0.2 


15* 


2.752 ± 0.040 


72.7 ±0.8 


K s 


17.3 ±0.5 


15= 


6.370 ±0.017 


209.4 ± 0.2 


.1 


13.8 ±0.4 


15 


6.352 ±0.017 


210.0 ±0.2 


H 


13.6 ±0.2 



a The position angle is measured from north to east. There is an additional 
systematic error of ~ 1 . 
^Apparent magnitudes were derived from TablefTl 

c The astrometry was measured relative to the photocenter of the inner binary. 
This object is listed in 2MASS with J = 1 1 .70 ± 0.06, H = 1 1.12 ± 0.08, and 
K = 1 1 .06 ± 0.06, and in DENIS with / = 12.43 ± 0.03. 

d This object is listed in 2MASS as having J = 14.48 ± 0.23, H = 13.81 ± 
0.32,andAT= 13.72±0.12. 

e The astrometry was measured relative to the south component B of the inner 
binary. 

^Marginal detection 

g Object located at edge of array. 

intensity by summing up all pixels within a circular aperture 
of 2" radius. 

The estimated Strehl ratios are shown in Figs.|3]&|4] except 
for the coronographic observations where we could not deter- 
mine ,R bs reliably, due to the presence of the mask. Because 
the coronographic targets are bright, and observed under sim- 
ilar conditions, we expect the Strehl ratios to be at the high 
end(>15%). 

The lower the Strehl ratio, the larger fraction of the stellar 
flux that is diluted into the seeing disc. Even down to Strehl 
ratios of a few percent, however, there usually is a diffraction- 
limited core of the PSF. This means that searches for point 
sources (such as stars) are greatly aided by AO even at low 
Strehl ratios, while searches for extended structures (such as 
circumstellar material) are critically dependent on high Strehl 
ratios. 

3.3. Contrast sensitivities 

In order to estimate the sensitivity to companions as a func- 
tion of separation, we used the following procedure: 

1 . For every pixel j in the detector, compute the distance 
rj to the determined center of the primary. 

2. Choose two radii Rq and Ri and fit an affine function 
f(r) = ciQ + a\r, where an, and a\ are fitting constants, to 
the intensity Ij of all pixels j with r, S [Rq,R\], Reject 
pixels more than 3a from the fit. 

3. Compute the standard deviation of the fit residuals, 
CpixCK) = StdDev[/j— f(rj)], where R is the mean of all 

n e [R ,Ril 

4. Measure the FWHM of the primary PSF, and the num- 
ber of pixels TVpsF and integrated flux F P s F within that 
area. 
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show the absolute sensitivity. The bottom horizontal axes show the angular separation and the upper axes the corresponding projected separation at the r\ Cha 
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lines depict the st ar/brown dwarf (0.0 8 M w : J = 12.5, H = 12.0, and K = 11.8) and brown dwarf/planetary mass (15Mj: J = 15.0, H = 14.6, and K = 14.2) 
boundaries for the Baraffe et al. 1 2003) evolutionary models at age 8 Myr and distance 97 pc. 
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VLT/ISAAC 




1.6 1.7 
Wavelength [um] 



FIG. 5. — ISAAC spectra of the primaries and companion candidates of 
ryCha 4, 9, 10, and 15, denoted in the figure. The spectra have been nor- 
malized and offset in steps of 0.25. The lowest spectrum shows the telluric 
atmospheric extinction, also normalized and of fset t o 0. All of these compan- 
ion candidates are likely background objects (§ 14.11 . 

5. The derived 5cr contrast sensitivity is then estimated to 
be 



£(/?) = -2.5 log 10 



5(7pix(fi)\/AfpSF 



•Fpsf 



mag. 



(1) 



The absolute sensitivity is obtained by adding the primary 
magnitude to £(/?). The function f(r) in step |2] is fit to re- 
move the radial gradient in the pixel intensities, that is due 
to the PSF. We confirmed that this procedure accurately esti- 
mates our achieved sensitivity by artificially placing intensity 
scaled PSFs at various separations. Since we limit the s earch 
to within 5", we do not correct for anisoplanatism (see § 13. 1> . 

For the semi-transparent coronographic observations, we 
made use of the measured suppression ratio in H (6.0 ± 
0.1 mag, §|2}. 

In the case of rj Cha 9, which is a tight 0."2 binary, we com- 
puted the contrast sensitivity from the center of light of the 
two stars, and estimated A/psf and fpsF as the average and 
sum, respectively, of both stars. The iso-intensity curves from 
the combined pair are elliptical rather than circular symmetric 
close to the stars, which is the reason the contrast sensitivity 
is estimated only at > 0."2 from the common center. 

The contrast sensitivities for all observations are presented 
in Figs.EI&IU 

3.4. Spectra of companion candidates 

The obtained ISAAC //-band spectra of the primaries and 
companion candidates of the systems 7/ Cha 4, 9, 10 and 15 




1.6 

Wavelength [urn] 

FIG. 6. — Resolved NACO spectra of the inner 0"2 r\ Cha 9 binary, normal- 
ized and offset by 0.25 from each other. The spectra are virtually identical. 
The scale has been chosen so that the flux scale / wavelength scale ratio is the 
same as in Fig. 151 

are shown in Fig. [5] Since the correction for telluric lines 
was poor due to high airmass, a plot of the extinction is also 
shown in Fig. [5] The NACO //-band spectrum of the inner 
0"2 binary of rj Cha 9 is shown in Fig.|6] 

4. DISCUSSION 
4.1. Companion candidates 

The -q Cha cluster is only 22° from the Galactic plane, so the 
probability of chance alignment stars in the 13"6x 13. "6 field 
of view is non-negligible. From the 2MASS All-Sky Data Re- 
lease we find that the density of stars brighter than H ~ 16 is 
~ 10~ 3 arcsec -2 , that statistically would produce about 3 ± 2 
chance alignments in our 17 fields, consistent with the total 
number of companion candidates (5) with H < 16 in our sam- 
ple. Since the observations are generally much more sensitive 
than this, we expect even more faint background stars. We 
have unfortunately no access to a deeper NIR survey of the 
region to find the local star density, but a general conclusion is 
that potentially all of the 10 new companion candidates found 
(Table|4j could be chance alignments. 

For four of the candidates (near 7/ Cha 4, 9, 10, and 6. "4 
from 77 Cha 15) we have NIR ISAAC //-band spectra. If these 
companion candidates had been members of the 77 Cha cluster, 
their age (^8 Myr) and luminosity (from Tableland the dis- 
tance 97 pc) would have implied v ery low-mass obje cts with 
atmospheres cooler than 2500 K (Baraffe et alj|200 3[). corre- 
sponding to spectral types later than M8.5 ( Kirkpatric k et alJ 
ll999t lLuhman et al. 2003). This is clearly not the case, as the 
ISAAC spectra reveal all companion candidates to have spec- 
tral types equal to, or earlier than, their primaries (Fig. [5jl. 
In particular, there is no evidence for water depression that 
should be visible at ~ 1.5 /im for spectral types later than M5 
dCushing et alJ20 05). We thus conclude that these companion 
candidates most likely are background stars, and not associ- 
ated with the T) Cha cluster. 

For the companion candidate near rj Cha 3, we have two 
epochs of data separated by 1 .25 yr that show a relative po- 
sition difference of (Aa,A6) = (48 ± 7,-14 ± 7) mas, im- 
plying a relative proper motion of (fi a ,[is) = (38 ± 6,-1 1 ± 
6) mas yr" 1 . A common proper motion is thus ruled out with 
> 5<7-significance, while the relative proper motion is roughly 
consistent with the companion candidate being a background 
star, since the proper motion of the 7/ Cha cluster is (/z a , jj,s) = 
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(-30.0 ±0.3, 27.8 ±0.3)masyr _I (Mamaie k~erall2000l) . 

The inner companion candidate of ?7Cha 3 also has two 
epochs of data, but unfortunately the positional precision is 
not sufficient to significantly constrain the relative proper mo- 
tion. The non-membership status of the star is instead re- 
vealed by its NIR colors; a H = 17.5 ?]Cha member should 
have had J-K ~ 3 rather than J-K = 0.1 ±0.7. 

The companion near candidate jyCha 1 is separated and 
bright enough to have 2MASS and DENIS photometry (Ta- 
bleg}. We reject this as a physical companion, since the ob- 
served I -(J ,H \K) colors ar e ~1 mag fainter than expected 
from the Baraffe et al. ( 1998) models, even when allowing for 
variability due to the epoch difference between 2MASS and 
DENIS. 

The remaining 3 companion candidate stars (near 7/Cha 
2, 7, and 8) have only single epoch //-band imaging, and 
therefore no direct way of ruling them out as physical com- 
panions. They are all quite distant from the system primary 
(> 5.6" ~ 540 AU) and faint (H > 16), however, and therefore 
likely chance alignments. 

We conclude that, among the 17 r/Cha members surveyed, 
there are no detected companions with projected separations 
20-500 AU. 

4.2. Orbit of riCha 1 

Kohl er & Petr-Gotzensl (120021) used multi-epoch observa- 
tions of r?Cha 1 between 1996-03-29 and 2001-12-10 to 
compute a preliminary orbit, estimating a dynamical mass. 
With the additional data point from 2004-04-03, we can con- 
strain the orbit further (see Fig. [7}. In particular, among the 
tw o preferred orbits with period s of 43 yr and 151 yr found 
in Kohler & Petr-Gotzens (2002), the longer orbit is clearly 
favoured by our data. To possibly find a better orbital solu- 
tion, we developed a simple orbit-fitting code that works as 
follows: 

1. For a given orbit k, compute the positions (sf,</>f) on 
sky for the dates of the observations, where s\ is the 
separation of orbit k and observation i, and e/>- the cor- 
responding position angle. 

2. Compute the square sum \ 2 of the differences between 
the positions (if,0f) from the assumed orbit, and the 
observed positions (s° hs , 0° bs ): 

where a s ,,■ is the error in separation and cr^ , the error in 
position angle (including estimated systematic errors). 

3. Find the orbit k that minimizes y^ik). 

The reason for using polar coordinates is that the position 
angle normally introduces a larger error than the separation, 
due to uncertainties in detecto r orientation. 

While the 151 yr orbit of Koh ler & Petr-Gotzensl (12001 
produces a good fit (Fig. [7}, we find a multitude of very dif- 
ferent orbits that produce equally good or better fits. In par- 
ticular, there is a general solution degeneracy such that ex- 
tremely eccentric orbits at high inclinations produce good fits 
(but with an unrealisticly high system mass). An example is 
given in Fig. [71 where the unbroken line shows an orbit with a 
period of 93 yr and semi-major axis 0"416, corresponding to 



0.2 




Aa farcsec] 

FIG. 7. — Relative astrometry between the two components of r) Cha 1. The 
crosses show the measured separations, with the right-most data point coming 
from this paper and the others from Kohler & Petr-Gotzens 12002). The Icr 
positional errors are indicated by the size of the cr osses. The clashed and 
dash-dotted curves are the orbit solutions found by Kohler & Petr-Gotzens 
1 2002, see their Table 4), with periods of 150.8 yr and 42.5 yr, respecti vely. 
The unbroken line is the unphysical 93 yr period orbit fit mentioned in § 14.21 
The small circles denote the positions in the orbits that correspond to the 
dates of the observations. 

the system mass of 7.6 M Q at the distance of 97 pc. Our con- 
clusion is that the orbit has to be followed for a longer time 
span before useful limits on the dynamical mass can be made. 
Alternatively, spatially resolved radial velocity measurements 
would add valuable constraints. 

4.3. Limits on companion probability 

The contrast sensitivity estimates from § 13.31 can be used 
to put limits on the number of likely companions. The basic 
approach is to assume that companions of stars are assigned 
by a stochastic process such that any given star system will 
have a companion with probability p, called the multiplicity 
(or binary) probability. The observed systems are then seen 
as a sample of this stochastic variable. Even if we knew with 
complete certainty that the 17 systems of ryCha had no wide 
stellar companions, there would still be a 5% chance that this 
outcome would have been produced with p = 0.16; thus, the 
95% confidence upper limit for p would have been 0.16. 

Before we can do the proper statistics, however, we need to 
correct for the observational biases. That is, given that there is 
a companion, would we have detected it? Let p be the proba- 
bility that a star has a companion, and qj the probability that a 
companion would have been detected in system j, given that 
there had indeed been a companion there. Then the probabil- 
ity that a companion is not detected in system j is l-pqj, and 
the probability that no companions are detected in N systems 
enumerated from 1 to yV is 

N 

®(p) = l[(l-pqj)- (2) 

7=1 

With a confidence set we mean the set of p such that 
$(p) < 1 ~~ a 7 where a is the chosen confidence, typically 
a = 95%. As <P(p) is a monotonously decreasing function of 
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p, the confidence set becomes an interval [0,p a ], where we 
call p a the a-confidence upper limit on p. 

To derive p a we need to know qj for all observed systems. 
For a given magnitude and separation from the star, this be- 
comes straightforward using the contrast sensitivity estimates 
from § 13.31 If the brightness is above the sensitivity limit at the 
given separation in system j, qj = 1; else qj = 0. Equation[2] 
thus reduces to $(/?) = (1 —p) n , where n is the number of sys- 
tems where detection would have been possible. This implies 
p a = 1 — (1 — a) 1 '". Because p a thus only depends on the dis- 
crete n, p a will also be discrete, as shown in Fig. [8] In the 
left panel of Fig. [8]we derive 95% upper limits on the prob- 
ability of companions of specific flux ratios to the primary, 
as a function of separation. To derive limits on the compan- 
ion probabili ty for spec i fic ma ss ratios q, we assume primary 
masse s from [Lyo et alJ (12 004). and use evolutionary models 
from lBaraffe et al.l J2003I) totranslate companion mass to an 
//-band magnitude, assuming the age 8 Myr, distance 97 pc 
and solar metallicity. In the right panel of Fig. [8] we instead 
display limits on the absolute sensitivity to companions of dif- 
ferent //-band magnitudes and masses. 

From a physical point of view, rather than knowing the 
companion probability as a function of observed (projected) 
separation, it is more interesting to find the probability as a 
function of semi-major axis. To do that we need to know the 
probability qj to find a companion, given that it has a specific 
semi-major axis. Our ability to detect a companion will de- 
pend on its projected distance to the primary, which in turn 
depends on viewing geometry and orbital phase (for eccentric 
orbits). In the Appendix, we compute the projected separation 
probability distribution for a companion of semi-major axis 
a. Assuming a random orientatio n, we derive the distribution 
analytically for circular orbits (§ IA.1> . and use a numerical 
approach in the case of eccentric orbits with an eccentricity 
distribution /(e) = 2e (motivated by both theory and obser- 
vations; see § IA.2> . The probability qj is then the probability 
that the companion is located at a detectable distance from the 
binary. Assuming that the sensitivity increases monotonously 
as a function of separation, this equals the probability that the 
companion is outside the "detection separation" d - i.e. 



1j = 



fs(s)ds=l-F s (d/d), 



(3) 



d/a 



where fs(s) is the projected separation probability de nsity dis- 
tribution found in the Appendix and shown in Fig. IA1I and 
Fs(s) is the corresponding probability distribution. Since the 
projected separation can be arbitrary small for any semi-major 
axis, we always have that Fs(d/a) > for d/a > 0, and thus 
qj < 1. Once we know qj, equation |2] is used to compute 
95% upper limits on the probability of companions. For the 
coronographic observations, the mask is at a fixed position 
only h = 2. "5 from the edge of the detector; we therefore intro- 
duce the additional correction factor g(s) = 1 - tt arccos(/i /s) 
(for s > h) into the integral of equation [3] where g(s) is the 
probability that a companion at separation s is in the field of 
view. 

In Fig.|9]we show the 95% upper limits in the same way as 
in Fig. [8] but as a function of the semi-major axes. 

4.4. Deficit of wide binaries 

It is remarkable that, despite our high sensitivity, the widest 
binary detected in the 77 Cha cluster is only 0."21, correspond- 
ing to the projected distance o f 20 AU. This lack of wide bi - 
naries was already noted bv iKohler & Petr-Gotzeiisl (|2002), 



albeit with smaller statistics and less sensitive measurements. 
Outside 20 AU, we would have found any stellar companion 
(of mass >0.07M©; see Fig. [8}. The fact that we do not de- 
tect any wide binaries among the 17 members implies that 
the wide (> 30 AU) binary probability w v cha with la con- 
fidence is lower than p\ a =0.01, and with 95% confidence 
is lower than /?95%=0.18 (Fig. [9}- This stands in stark con- 
trast to the other nearby, young (of the similar age 8 Myr) 
TW Hydrae association (TWA), where 11 out of the mem- 
ber stars TWA 1-19 a re binaries with separations > 30 AU 
jBrandeker et alJ2003l) . corresponding to a wide binary prob- 
ability of wtwa = 0.58^o [4, clearly different from what we find 
in 77 Cha. This difference cannot be due to detection sensitiv- 
ity differences; all the companions in TWA would have easily 
been detected in the 77 Cha cluster by this survey. The differ- 
ence cannot be explained by small number statistics either, as 
the statistical likelihood that the 77 Cha cluster and TWA have 
binary frequencies drawn from the same distribution is less 
than 2 x 10~ 4 (see the Appendix). 

There is certainly no deficiency of close binaries (< 30 AU) 
in the 7/ Cha cluster dLvo et alJ2004l) . If anything, there might 
be a small over-abundance, compared to TWA, but the statis- 
tics are inconclusive. 

One may speculate on the cause for such a difference in 
the wide binary population between the 77 Cha cluster and the 
TWA. One notable difference is that the 77 Cha cluster is much 
denser than the TWA. Could it be that wide binaries are dy- 
namically disrupted in the 77 Cha cluster? To investigate this 
possibility we estimate the timescale for a binary to undergo 
a strong encounter with another star in the cluster. From 
Ivano vaet alJ (12005, equation 7), we have that this timescale 



r coll = 1.7 x 10 8 yr x 77 2 /T 2 «;' ^ 



M\M\ 



2 M\ +Mo+(M) , , 

x I+77 — — (Af) 

1 k MiM 2 



(4) 



where 77 is the hardness of a binary, k ~ 2, 715 is the number 
density of star systems in units of 10 5 pc~ 3 , (Af) is the mean 
mass per star system in units of M Q , and M\ and M2 are the 
masses of the binary components, also in units of M Q . The 
hardness of a binary is defined as 



MiAfz 

\M)o l a 



(5) 



where G is the gravitational constant, a is the binary sepa- 
ration, and a is the velocity dispersion of the cluster. Bina- 
ries with 77 < 1 are termed soft, and those with 77 > 1 hard. 
In general, soft binaries are disrupted by strong encounters, 
while hard binaries may survive. In the 77 Cha cluster, we 
have the cumulative m ass 16.6 M Q distributed over 17 sys- 
tems (lT,vo et al 12 004). giving (A/) = 1 .0. There are presently 
no accurate radial velocities published for the 77 Cha mem- 
bers. However, from the estimated age of the cluster (~8 Myr) 
and observed effective radius (~ 0.2 pc), we infer that the 
cluster most likely is gravitationally bound. The present 
escape velocity is namely v esc = (2GM/R) X I 2 ~ 0.75 km s" 1 , 
assuming M = 13M Q within R = 0.2 pc dLvo et al J 12004). 
which gives a crossing time of merely 0.3 Myr, enough to 
traverse the cluster core 25 times during its lifetime, while 
the relaxation time is on ly a few times the crossing time 
( iBinnev & Tremainelll987L chapter 4). We therefore assume 
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FIG. 8. — The 95% confidence upper limits on the probability of companions. The left panel shows the upper limits for companions of specific contrasts to the 
primary, as a function of separation from the star. The dashed lines (colored blue in electronic edition) show the sensitivity limits to mass ratios q = M com p/M pl j m , 
while the unbroken lines (colored red) show the limits for flux-ratios AH. The right panel displays the absolute sensitivity, where the dashed lines (colored blue) 
show sensitivity to companion masses, and the unbroken lines (colored red) show the sensitivity to apparent //-band magnitudes. The dashed horizontal line 
shows the 95% confidence upper limit that would have been obtained statistically for our sample of 17 systems, had the detection sensitivity been 100%. The 
increased upper limits at separations >2"5 is due to part of the field being outside the array. 
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FIG. 9. — Same as in Fig. |8| except that the sensitivity limits are displayed as a function of semi-major axes, assuming the eccentricity probability density 
distribution f(e) = 2e. 



o — Vesc/2 = 0.37 km s" 1 , in accordance with the virial the- 
orem, implying that the condition for a binary to be soft 
in the ?7Cha cluster is a > 1600AU > 30 AU for the typi- 
cal masses Mi = M2 = 0.5. Moreover, for the r/Cha cluster 
«5 ~ 3 x 10~ 3 , and the collision timescale for 77 = 1 binaries is 
thus t co ii ~ 26Gyr, i.e. 3000 times longer than the lifetime of 
the system. That dynamical interactions between binaries in 
the present configuration of the cluster should be responsible 
for the lack of wide (>30AU) binaries, consequently seems 
highly unlikely. 

Since both the 77 Cha cluster and TWA are of similar age, 
the remaining explanation is that the difference in multiplic- 
ity properties were imprinted during the formation phase, as a 
result of different initial conditions. Either the groups formed 
with different multiplicity properties, or the properties dy- 
namically evolved very early on, when the stars were possibly 
much closer together. 

There seems to be a general trend that denser groups 



have smaller wide binary frequencies than sparser regions; 
the sparse regi ons Taurus, Ophiuch us and Chamael eon 
JDuchend 11991 . and MBM 12 JBrandekeret all 12001 all 
have high wide binary frequenci es, while the denser re- 
gions Trap ezium JPetr et alJll998l) . and NGC 2024, 2068, 
and 2071 (Padgett et al. 1997), have low wide binary fre- 
quencies. Solar-type main-sequence stars in th e solar neigh- 
borh ood have a binary probability of ~0.45 ( iLeinert et alJ 
1 19931 and references the rein), with half being wide 
(Dua uennov & Mavoill 19911) . and thus wms = 0.23, which is 
right in between w r) cha and wtwa- It is therefore not clear 
from the wide binary statistics alone if sparse or dense star 
formation is the dominant mode - possibly both contribute 
equally. What is clear is that models of cluster formation and 
early evolution probably are essential to exp lain multiplicit 
pro 

iCoodwin et all2006). 



properties (e.g. Bate et al. 2002; Delgado-Donate et al. 200- 



z 



5. CONCLUSIONS 
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We summarize our conclusions as follows: 

1 . We found no new companions to stars in the rj Cha clus- 
ter, despite being sensitive down to the star/BD limit 
outside 0."3 (30 AU) and down to the BD/planet limit 
outside 0."5 (50 AU). 

2. We have constrained the orbit of ryChal further, but 
are unable to usefully constrain a dynamical mass. Re- 
solved radial velocities of the two components, or a 
longer astrometric time baseline, are required for an ac- 
curate mass estimate. 

3. The 95% upper limit for the wide (>30AU) binary 
probability w,,cha in r?Cha is /?95%=0.18. This con- 
trasts to the wide binary probability Wtwa = 

of TWA. The likelihood that w^cha = wtwa * s ' ess tnan 
2 x 10" 4 . 



4. Multiplicity properties depend on the initial conditions 
of the formation environment. 
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APPENDIX 

A. PROBABILITY DISTRIBUTIONS FOR PROJECTED BINARY SEPARATIONS 

This appendix concerns the computation of probability distributions of the observed projected separation between the two 
components of a binary system, given that the semi-major axis is known. Here we use the semi-major axis a of the companion 
orbit relative to the primary. If the mass of the companion is a significant fraction of the primary, this will be different from 
the semi-major axis a com of the companion orbit relative to the center of mass. The relation is a = (M comp /M pr i m + l)a CO m, where 
M comp and M pr ; m are the masses of the companion and primary, respectively. The derived probability distributions below are 
easily scaled to a com in case of specific mass-ratio systems. 



A.l. Companions in circular orbits 

To compute the probability distribution for the projected separation of a companion, we need to make the assumption that the 
orbit is observed from a uniformly distributed direction; that is, any viewing direction is equally probable. 

To generate a stochastic vector with a direction uniformly distributed over the unit sphere is a less trivial problem than one 
might naively expect. One cannot simply use a uniform distribution of spherical coordinates, since this would bias the vectors 
too much towards the poles; and one cannot simply use uniformly distributed Cartesian coordinates, since this would bias the 
directions towards the corners. Instead, a stochastic vector d = (x,y,z) uniformly distributed over the unit sphere may be generated 
as 

x= \J 1 -z 2 cos^ 

y= \/ 1 -z 2 sin</!> 
zeU(-l,l), 

where X £ U(a,b) means that X is a stochastic variable uniformly distributed between a and b, and 4> G U(0,2w). This is a 
corollary from a theorem by Archimedes, that states that a lateral area of a section cut out of a sphere by two parallel planes 
equals A = 2nRh, where R is the radius of the sphere and h is the distance between the planes. 

With R being the radius vector from the star to the companion, R its length and d the unit viewing direction vector, the projected 
distance between the star and the planet onto a plane perpendicular to the viewing direction will be s = ^R 2 -(K»d) 2 . 

Since d is uniformly distributed on the unit sphere, we can (without loss of generality) let R = Re z , where e z is the unit vector 
along the z-axis. The scalar product becomes R»d = Rd z , where d z is the z-component of the direction vector d. That is, the 
scalar product is only dependent on the z-component of the direction vector: 

s = R^\-(%. (Al) 

The probability distribution for the corresponding stochastic variable S = V 1 — D 2 , where D 6 U(0, 1 ), becomes Fs(s) = 1 - vl -s 2 , 
and the probability density distribution fs(s) = dFs(s)/ds = s/y/l-s 2 . 

These distributi ons t hus describe the projected separation s, in units of the orbital radius, of a companion in a circular orbit 
around a star (Fig. lAU . 



A.2. Companions in elliptical orbits 

To derive the projected separation distribution for elliptical orbits is slightly more complicated than for circular motion, since 
the radial distance between the star and its companion is non-linear in time. We also have to make an additional assumption on 
the eccentricity distribution. We still do not have to bother with the orbital elements concerning a specific orientation of an orbit 
however, since we already have assumed that the viewing direction is random. 
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FIG. Al . — The left panel shows the probability density distribution for the projected separation of a companion to a star, in units of its semi-major axis. Two 
cases are plotted, one where the orbits are assumed to be circular (solid line), and one where the eccentricity density distribution is assumed to be f(e) = 2e 
(crosses). The projected separatio n distr ibution is computed exactly for circular orbits, and by a monte carlo approach for the eccentric orbits. The dashed line is 
an approximate fit, as outlined in § IA.2I The right panel shows the corresponding probability distribution, i.e. the integral over the probability density. 



The mean anomaly M of an elliptical orbit is defined to be linearly increasing with time <Meeuslll9911 Chapter 29). Its 
relation to the eccentric anomaly E, that is used to compute the actual position of the companion, is defined by the Kepler 
equation E = M+esinE, where e is the eccentricity of the orbit. Since the Kepler equation is transcendental it cannot be solved 
algebraically, which complicates the analysis. We therefore derive the probability density distribution using the following monte 
carlo approach instead. 

1. Let the mean anomaly M 6 U(0,2ir). 



2. Get the eccentricity e from some pre-defined distribution /(e). 

3. Solve Kepler's equation numerically and compute the eccentric anomaly E(M). 

4. Compute the instantaneous distance between star and companion, R = 1 -ecosE. 



5. Let the projection term d z £ U(-l, 1), in line with equation lAll and compute the projected separation s = R\/\ ~d\. 



The choice of ec centricity distribution is important for the resulting probability density distribution. Theoretical considerations 
( Ambartsu rmanlfl937l) predict the distribution f(e) = 2e, where e € [0,1], which is reasonably confirmed by observations of 
long-period binaries (>1000 days; Duquennov & Mavoi 1991). To generate a st ocha stic variable with f(e) = 2e, we let e = X l l 2 , 
where X £ U(0, 1). The result from a simulation of 10 8 binaries is shown in Fig. lAll The probability density distribution is well 
approximated, although not perfectly, by the function 



7T /7T \ 

f s (s) = -sm^-s) , 

where s £ [0,2], which is over plotted in Fig. lAll The corresponding probability distribution is 



'7T 

-cos I — s 



(A2) 



(A3) 



A stochastic variable S with this distribution can be generated from I G f(0, 1) by setting S = 2ir 1 arccos(l - 2X), that might 
prove useful for future completeness studies. 



B. MULTIPLICITY STATISTICS 

Statistical estimates of multiplicity frequencies are often limited by small sample numbers, where approximate methods using 
assumptions of Poisson or normal statistics become insufficient. Instead, more accurate estimates can be obtained by explicit 
use of binomial statistics. The multiplicities of N systems of a stellar association are viewed as outcomes Xi of a binomially 
distributed stochastic variable X £ Bin(l,p), where p is the multiplicity probability for any specific system, x,- = 1 if system i 
is multiple, and x, = if not. The number of multiple systems is then k = Yli x i< which m itself can be seen as the outcome of 
a binomially distributed stochastic variable K £ Bin(N,p). The goal is to constrain p from a given set {x,}, assuming all x,- are 
drawn from the same distribution. Here we summarize results from two outstanding problems, how to constrain the multiplicity 
probability p by a confidence interval, and how to decide if two sets {x,} and {y,} are outcomes of the same binomial distribution. 
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B.l. Binomial confidence intervals 

The well-known maximum-likelihood estimator of p is p = ^tXi/N, but to compute an accurate variance of this estimator is 
much more difficult. The most widely used way of computing the standard deviation of p is by using the so-called Wald method, 
that has the simple form a = (p(l - p)/N) x l 2 . The Wald standard devia tion produces ver y poor approximations whenever p is 
close to or 1, however, and its usage is generally not recommended (B rown et alJl200lT) . The problem of computing accurate 
confidence interva ls for the binomial distribution is not a new one, and there exist a plenitude of literature on the subject (see 
iBrown et alJ200l| for a review of methods). A robust "exact" method guaranteed to produce intervals with confidence of at least a 
was proposed bv lClopper & Pearsonl J1934 . If K € Bin(N, p), the probability that an outcome is k is P(k of N)= ( f [)p k (l- p) N ~ k , 

and the probability that k < n < N is P(k < n of AO = Y^k-o^^ °f F° r an y g lVen significance a and specific outcome k, the 
Clopper- Pearson method consists of finding the interval (/?mi n ,Pmax) such that if p > p^, then P(j < k of AO < (1 - a)/ 2, and if 
P < Pmax, then P{j > k of N) < (1 - a)/2. By solving for the equalities, we get 

i:( N )pL^-p^f- j ~ cbi) 

1=0 KJ/ 

E l fN\ : „ ; l + a 

. )PU(1-P^f- J = ^, (B2) 

where we have used that P(j > k of N)= l—P(j < k+ 1 of AO- The equations IB 1 1 & lB2l are normally best solve d num erically. 2 
In the degenerate case where k = or k = N, the confidence interval becomes one-sided, and the equations IB 1 1 & lB2l are easily 
solved analytically (e.g., equation|2]in § 14.31 1. 

B.2. Comparing binomial distributions 

Let &a out of Na systems in association A be multiple, and ks out of Nb systems in association B. The question that naturally 
arises is, is the multiplicity probability pa of system A similar to pr of sy stem Bl As in the case of binomial confidence intervals, 
there are several tests available in the literature (see Storer & Kim 1990 for a review). A conservative "exact" hypothesis test is 
based on Fisher ( 1935), where the hypothesis that pa = pb is tested: Let h(kA,kB,NA,NB) = ( k f) / (j^+fo) an d I[expression] be 
the indicator function that is 1 if expression is true, and otherwise. Then the test function is 

mm(n A ,k A +k B ) 

T= 22 h(x,k A +kB-x,nA,n B )xI[h(x,k A +kB-x,nA,nB) <KkA,k B ,n A ,nB)], (B3) 

x=max(0,kA+kB— /is) 

and the hypothesis is rejected if T < a, where a is the significance of the test. As an example, if k A = 0, Na = 17, k B = 1 1 , and 
Nb = 19, then T = 1 .46 x 10~ 4 , which is the quoted likelihood that the wide binary probability is equal in r/ Cha and TWA (§ I4.4I >. 
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